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Abstract: We develop a new approximation scheme aiming at extracting higher-point corre- 
lation functions from the JIMWLK evolution, in the limit where the number of colors is large. 
Namely, we show that by exploiting the structure of the 'virtual' terms in the Balitsky- JIMWLK 
equations, one can derive functional relations expressing arbitrary n-point functions of the Wil- 
son lines in terms of the 2-point function (the scattering amplitude for a color dipole). These 
approximations are correct not only in the regime of strong scattering, where the evolution is 
indeed controlled by the 'virtual' terms, but also in the regime of weak scattering, where they 
reduce to the corresponding BFKL solutions. This last feature follows from the fact that the 
JIMWLK Hamiltonian is a linear combination of the pieces responsible for the 'real' and 'virtual' 
terms, respectively. We apply this scheme to two examples: the 'color quadrupole' (the 4-point 
function of the Wilson lines which enters the cross-section for the production of a pair of jets 
at forward rapidities) and the 'color sextupole' (the 6-point function). For particular configura- 
tions of the quadrupole, our general formula reduces to relatively simple expressions that have 
been previously proposed on the basis of the McLerran-Venugopalan model and which were 
recently shown to agree quite well with exact, numerical, solutions to the JIMWLK equation. 
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1 Introduction 

Since the derivation, more than a decade ago, of the equations describing the high-energy evo- 
lution in QCD to leading logarithmic accuracy — the infinite hierarchy of Balitsky equations 
[1] for the n-point functions of the Wilson lines (= scattering amplitudes in the eikonal approx- 
imation) and the functional JIMWLK 1 equation [2-9] for the color glass condensate (CGC = 
the small-x part of the wavefunction of an energetic hadron) — , there has been little progress 
towards solving these equations for correlations which are more complicated than the 2-point 
function (the scattering amplitude of a color dipole). The progress has been mostly hindered by 
the extreme complexity of these equations (an issue that we shall shortly return to), but also 
by the fact that, for quite some time, the dipole amplitude was the only such a correlation to 
be directly relevant to phenomenology (e.g., for deep inelastic scattering and for single-inclusive 
particle production in hadron-hadron collisions; see e.g. [10-12] and references therein). 

However, the situation has changed in the recent years, with the advent of new, less inclu- 
sive, data, which probe higher-point functions via the multiparticle correlations in the particle 
production at high energy. In particular, the recent data at RHIC, on the azimuthal correlations 
in the forward di-hadron production in deuteron-gold collisions [13], are sensitive to a 4-point 
function of the Wilson lines known as the 'color quadrupole' [14-17]. These data, together with 
their successful description by a phenomenological analysis using the general ideas of the CGC 
[18], have revived the interest in the higher-point correlations and triggered new efforts in that 
sense, mostly in relation with the quadrupole [19-22]. But these efforts have been hampered 
by the main difficulty alluded to above — the extreme complexity of the Balitsky- JIMWLK 
equations — , which so far has precluded the obtention of any analytic solution, even approxi- 
mate. It is our purpose in this paper to present a major step in that sense, in the form of an 
approximation scheme for the higher-point functions with n > 4, that we believe to be new. 
But before describing our approach, let us briefly recall the state of the art. 

1 'JIMWLK' stands for Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov, Kovner. 
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The only analytic estimate for the quadrupole which is currently available within the CGC 
framework is that obtained within the McLerran-Venugopalan (MV) model [23, 24], which refers 
to a large nucleus at not too high energy. This estimate has first been computed in the limit 
of a large number of colors N c in Ref. [14], and more recently for generic N c in Ref. [20]. The 
evolution equation for the quadrupole is also known: this has first been derived in Ref. [14] at 
large N c , by using the dipole picture [25], and then by several authors [21, 26, 27] including 
ourselves, for generic N c , using the JIMWLK Hamiltonian. (Our version of the derivation is 
presented in Appendix A.) But this is only one out of an infinite hierarchy of coupled equations, 
which moreover exhibit a complicated non-local structure in the transverse coordinates. This 
non-locality becomes more and more severe with increasing the number n of external points. 

The Balitsky- JIMWLK hierarchy considerably simplifies in the large- N c limit, where it 
reduces to a triangular hierarchy. It is then enough to solve a finite number of equations in 
order to determine a given correlation function. The dipole obeys a closed, non-linear, equation, 
known as the Balitsky-Kovchegov (BK) equation [1, 28]. The quadrupole evolves according to 
an inhomogeneous equation, in which the dipole plays the role of a source term. The equation 
for the sextupole (a color trace of six Wilson lines) also involves the dipole and the quadrupole, 
and so on. But even in that limit, the non-locality of the equations is such that it renders 
prohibitive their (numerical or analytic) study, except in the simplest case: the BK equation, 
which has been studied at length [10-12]. As a matter of facts, the most promising approach 
towards numerical studies refers to the functional JIMWLK equation by itself, and not to the 
ordinary evolution equations in the hierarchy: this functional equation can be reformulated as 
a Langevin process [29], which is well suited for the numerics. The feasibility of such numerical 
solutions has been demonstrated in Refs. [30-32], which however focused on the 2-point function 
(the dipole amplitude) alone. It was only very recently — when the present study was under 
progress — that this numerical method has been extended to the quadrupole [22]. 

The analysis in Ref. [22] follows the evolution of the quadrupole with increasing energy, 
starting with initial conditions of the MV type and for two special configurations of the four 
external points in the transverse plane: the 'line' and the 'square'. (These configurations are 
illustrated in Figs, l.a and b below.) Interestingly, the results show a good agreement with a 
heuristic extrapolation to high energy of a formula, in particular an expression for the quadrupole 
in terms of the dipole, that was derived within the MV model [14, 20]. Such an agreement is hard 
to explain without a more fundamental understanding of the relation between the extrapolation 
allude to above and the actual JIMWLK dynamics. Moreover, the continuation of the numerical 
solutions in [22] to more systematic studies (say, in view of the phenomenology) can be very 
difficult and demanding in computer power. Thus, it is necessary to remedy such shortcomings 
by completing the numerical solutions with controlled analytic approximations — a task that 
we shall accomplish here. 

Specifically, our strategy will consist in solving a simplified version of the Balitsky-JIMWLK 
equations at large N c , in which we keep the 'virtual' terms, but we drop the 'real' ones. The 
distinction between 'real' and 'virtual' is meant here in the same sense as in the dipole picture 
[25], that is, it refers to the evolution of the projectile (dipole, quadrupole, etc): the 'real' terms 
describe processes in which the projectile splits via the emission of an additional, soft, gluon 
(at large N c , this yields e.g. a color dipole splitting into two dipoles), whereas the 'virtual' 
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terms refer to the complementary probability that the projectile survive without splitting. By 
themselves, the 'virtual' terms control the evolution of the projectile S-matrix in the vicinity 
of the unitarity (or 'black disk') limit, so our approximations are a priori justified in the strong 
scattering regime deeply at saturation. But our scope in this paper is more general than that: 
by exploiting the structure of the 'virtual' terms, we shall derive explicit expressions for the 
quadrupole and the sextupole in terms of the dipole, which represent global approximations, 
valid for both weak and strong scattering. We shall indeed verify that, in the regime where 
the scattering is weak, our general results reduce to the expected, linear, relations between the 
quadrupole, or the sextupole, and the dipole scattering amplitude. This is not an accident: it 
follows from the fact that a linear relation is always preserved by an evolution equation derived 
from a Hamiltonian — in particular the 'virtual' part of the JIMWLK Hamiltonian. (We recall 
here that the JIMWLK Hamiltonian is the direct sum of the pieces responsible for the 'real' and 
'virtual' terms, respectively.) Hence, by using the BFKL approximation for the dipole scattering 
amplitude within our general expression, we recover the correct BFKL results for the quadrupole 
and the sextupole (although the BFKL evolution is sensitive to both 'real' and 'virtual' terms, 
of course). 

Our method is general and systematic: by using similar techniques, it is possible to derive 
expressions for all the n-point functions of the Wilson lines in terms of the dipole <S-matrix. 
Furthermore, our general expressions are also consistent with initial conditions of the McLerran- 
Venugopalan type (at large N c ). Hence, they provide a unified description of the initial conditions 
and of their high-energy evolution, in both the dilute (BFKL) and dense (saturation) regime. 
In our opinion, the ultimate reason why such a strategy can work is the fact that, within the 
JIMWLK evolution, the multi-gluon correlations are built exclusively via high-density effects, 
that is, via gluon recombination in the approach towards saturation. 

The numerical analysis in Ref. [22] provides already a test of our approximation scheme: 
for the particular quadrupole configurations investigated there (the 'line' and the 'square'), 
our general result coincides with the high-energy extrapolation of the respective MV formula?, 
that in Ref. [22] were found to agree quite well with the numerical solutions to the JIMWLK 
equation. The agreement found there is even better when using the finite-A^ version of the 
MV formula?, showing that the finite- N c corrections to the evolution of the quadrupole yield 
somewhat sizeable effects when N c = 3. This should be contrasted with the corresponding 
situation for the dipole, where one has numerically found [30, 31] that the finite- N c corrections 
in the JIMWLK evolution are surprisingly small — at the level of the percent accuracy instead 
of the 10% (~ 1/Nc with N c = 3) that would be normally expected. 

Let us finally mention some possible limitations of our present analysis, which could be 
improved by further studies. A priori, our approximations are better justified for relatively 
symmetric configurations of the n-point functions, which are such that the transverse separations 
between the opposite-sign charges (i.e. between quarks and antiquarks) are not very different 
from each other. It would be interesting to test this limitation by comparing our predictions for 
very asymmetric, 'small-large', configurations — which turn out to be very interesting (see the 
discussion in Sect. 4) — to numerical solutions to the JIMWLK equation. We also note that, by 
construction, our approximations are guaranteed to be correct in the weak scattering regime and 
in the approach towards the black disk limit, but not necessarily also in the transition between 
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these two regimes, as occurring around the saturation scale. Still, our results provide a smooth 
(infinitely differentiable) interpolation between the two limiting regimes and as such we believe 
them to be qualitatively and even quantitatively correct including in the transition region. This 
is again in agreement with the numerical analysis in Ref. [22]. 

The paper is organized as follows. In Sect. 2 we present the evolution equations for the 
dipole and the quadrupole that we shall use as examples to illustrate general properties of 
the Balitsky-JIMWLK hierarchy, like the relation between 'real' and 'virtual' terms and the 
simplifications which appear in the large- N c limit. (The respective equations will be explicitly 
derived in Appendix A and this will give us the opportunity to describe the method that we 
shall later use, in Appendix B, to construct the corresponding equation for the sextupole.) In 
Sect. 3 we proceed with a study of the evolution deeply as saturation, as described by the 
virtual terms, and introduce our approximation scheme on a simple, pedagogical, example: the 
'line' configuration for a quadrupole. In Sect. 4 we construct our global approximation for a 
generic configuration of the quadrupole and study some limiting cases corresponding to simple, 
but interesting, configurations. In particular, we shall identify configurations for which the 
quadrupole is exactly factorizing in the product of two dipoles. The corresponding analysis of 
the sextupole, including its evolution equation, is given in Appendix B. 



2 Evolution equations 

In this section, we shall review the Balitsky-JIMWLK evolution equations for the dipole and 
the quadrupole and explain the origin and the physical significance of the various terms which 
appear in these equations. This discussion will be useful in view of the construction of an 
approximation scheme later on. 

Within the CGC effective theory, valid to leading logarithmic accuracy at high energy, the 
evolution of a gauge invariant operator O with increasing energy is determined by 

d -§ L = (HO)y, (2.1) 

where the brackets refer to the average over the color fields in the target, as computed with the 
CGC weight function, and H is the JIMWLK Hamiltonian. The latter, when acting on gauge 
invariant observables, is most conveniently given in [33] and reads 

H ~ "is? L + ^ " n% ~ m Y <2 2) 



In the above equations, a a {x ,x) is the target color gauge field in the covariant gauge, where 
Aa = 5 M+ a a , M is the dipole kernel 

M = (tt u)^ 

and and V are Wilson lines in the adjoint representation: 

Vl = Pexp (ig I &x-a a {x-,x)T a \ , (2.4) 



with P denoting path-ordering in x . Our conventions are such that the nuclear target (the 
CGC) is a right-mover, whereas the projectile (to which refers the operator O) is a left-mover. 
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The functional derivative in Eq. (2.2) are understood to act at the largest value of x~ , that is, 
at the upper end point of the path-ordered exponential in Eq. (2.4). 

The operators that we shall deal with are all built with products of Wilson lines (in the fun- 
damental representation, for definiteness) and represent the 5-matrices for systems of partons 
(quark and antiquarks) scattering off the nuclear target, in the eikonal approximation. Specif- 
ically, we shall focus on the color dipole — a quark-antiquark pair in an overall color singlet 
state, with S-matrix operator 

S X ix% = SxiX2 = j-j- tr(V r a [ 1 K K2 ) , (2-5) 

the color quadrupole — a system of two quarks and two antiquarks, for which 

Qx!x 2 x3X4 = S X1X2 X3X4 = ^y - tr(K 1 f 1 VJc 2 V^ 3 V^ 4 ) , (2-6) 



and the color sextupole, with 



Sxix 2 x s x 4 x 5 x 6 — tr(V^ V^ 2 V^ 3 V^V^VJcg). (2-7) 

Higher-point correlators can be similarly considered, including those which involve the product 
of several color traces (see e.g. Eq. (2.12) below). In Eqs. (2.5)— (2.7), and V are Wilson lines 
in the fundamental representation. The action of the functional derivative on these Wilson lines 
reads 

Kl = igS xu t a Vl -^-V x = -igSxuVx t a . (2.8) 



Sa a x x 5a a 

By using these rules within Eqs. (2.1)— (2.2), it is straightforward to derive the evolution equa- 
tions satisfied by the expectation values of the three operators introduced above. A streamlined 
derivation of the respective equations for the dipole and the quadrupole will presented in Ap- 
pendix A, while the corresponding equation for the sextupole will be shown (without the details 
of the derivation) in Appendix B. 

The ensuing equation for the dipole looks relatively simple (we denote a = a s N c /ir): 

d{S Xl X 2 }Y _ 0_ f KA /A q _ q \ /r, Q\ 

„ / ■> V lXlX2Z\>- , XlZ>-'ZX2 '-'X1X2/Y 1 V^' i '/ 



dY 2vr . 

but that for the quadrupole is considerably more involved: 



d(Qx 1 X2X S X 4 )Y _ 0_ f I", . A _i_ \A KA \lq A \ 

QY J \- /V lxiX2Z ~T J V *-X\X4Z J Vl X2X4Z) \)JX1ZW ZX2X3X4 / Y 

{A // txiX2Z A4 X2X3Z J^X\XiZ ) (S Z x 2 QxiZX$X4 /Y 

+ (M X2X3Z + M X3X4Z -M 

X2 X4Z ) {SxszQxi X2ZX4) Y 
+ (M X4X4Z + M X3X4Z -M X1X3Z ){^ZX4 Qx\X2X%z)Y 

{•^ // ^xiX2Z A^ajg^z + ■A / ixjX4Z ~i~ ■^ / ^X2X3z){Qxix 2 X3X4)Y 
(A / J a ; 1 a; 2Z + ■A / txsx 4 z J^xiXiZ J^X2X4z) {S*xiX2 SxsX4lY 
- (M X4X4Z + M X2X3Z -M xix-jz -M X2X4z){Sx 3 X2^X 1 X4)Y -(2.10) 
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We shall now discuss the role of the various terms in the above equations, in order to explain, 
in particular, the difference between 'real' and 'virtual' terms. 

Consider first Eq. (2.9) for the dipole S'-matrix. The 'real' term is the first term in the right 
hand side, which is quadratic in S. When interpreted in terms of projectile evolution, this term 
describes the splitting of the original dipole (a?i, X2) into two new dipoles (a?i, z) and (z, X2), 
which then scatter off the target. The 'virtual' term, that is, the negative term which is linear in 
S, describes the reduction in the probability that the dipole survive in its original state — that 
is, the probability for the dipole not to split. The relative size of these two types of terms ('real' 
and 'virtual') is constrained by probability conservation, as correctly encoded in the JIMWLK 
Hamiltonian. Remarkably, the 'real' term has been fully generated by the last two terms in 
the Hamiltonian (2.2), whereas the 'virtual' term originated from the first two terms there. 
All the terms appearing in Eq. (2.9) are of leading order in N c . Subleading terms, of relative 
order 1/N%, had been also generated, at intermediate steps, by each of the four terms in the 
Hamiltonian, but they have exactly canceled after summing up all the contributions. Note the 
cancelation of 'ultraviolet' (i.e. short-distance) singularities between the 'real' and the 'virtual' 
terms: the dipole kernel (2.3) has poles at z = u and z = v, which within the integral over z 
generate logarithmic divergences, separately for the 'real' and for the 'virtual' term. However, 
these divergences cancel in the overall equation, because of probability conservation together 
with the property of 'color transparency' : S Xl x 2 ~* 1 when x\ — > xi- The latter is merely the 
statement that a 'color dipole' of zero transverse size is a totally neutral object, which cannot 
interact. 

An entirely similar discussion applies to the evolution equation (2.10) for the quadrupole. 
The terms involving (SQ)y in the right hand side are 'real' terms describing the splitting of 
the original quadrupole into a new quadrupole plus a dipole, and have been all generated by 
the action of the last two terms in the Hamiltonian (2.2). The 'virtual' terms involving (Q)y 
and (SS)y are necessary for probability conservation, and have been generated by the first two 
terms in the Hamiltonian. Once again, all the terms subleading at large N c (as generated by 
the individual pieces of the Hamiltonian) have canceled in the final equation. Furthermore, all 
possible ultraviolet divergencies due to the poles of the dipole kernel cancel between the 'real' 
and the 'virtual' terms, due to color transparency. 

All the above features are generic: they also apply to the Balitsky-JIMWLK evolution of 
the sextupole and, more generally, to any (gauge-invariant) operator which involves a single 
trace of products of pairs of Wilson lines, of the form 

S Xl X2...X2n-lX2n = ^ T ( Vsei ^2 ■ • ■ ^£D2«-1 ^ x 2n ) • (2-11) 

As for the multi-trace operators, of the form 




(2.12) 



the only new feature is that the respective evolution equations involve genuine corrections, 
as generated when the two functional derivatives in Eq. (2.2) act on Wilson lines which belong 
to different traces (see e.g. Appendix F in [34] for an example). 

As manifest by inspection of Eqs. (2.9) and (2.10), these equations are generally not closed: 
e.g., the equation for the quadrupole also involves the 4-point function (SS)y and the 6-point 



- 6 - 



function {SQ)y, which in turn are coupled (via the respective evolution equations) to even 
higher-point correlators, so altogether one has to deal with an infinite hierarchy of equations. 
But important simplifications occur in the large-iV c limit, as anticipated in the Introduction. 
Then, for a multi-trace operator like (2.12), it is not hard to show that the hierarchy admits 
the factorized solution 

(O)y = (i^ te (^i^--0> K (^te(^^«--0> K (^M^4% ! i---)) v ---, (2-13) 

provided this factorization is already satisfied by the initial conditions. Then the hierarchy 
becomes triangular and therefore finite: Eq. (2.9) becomes a closed equation for (S)y (the BK 
equation), while Eq. (2.10) becomes an inhomogeneous equation for (Q)y with coefficients which 
depend upon (S)y- This is the equation originally derived in [14] and that we shall further study 
in the next section. 



3 Prom virtual terms to global approximations 

In order to better appreciate the approximation scheme that we shall eventually propose, it is 
useful to investigate the behavior of the solutions to the evolution equations introduced in the 
previous section in two limiting regimes: the dilute, or weak scattering, regime, where all the 
transverse separations r%j = \xi — Xj\ between the external points are very small compared to 
1/Q S (Y), and the dense, or strong scattering, regime, where all these distances are much larger 
than 1/Q S (Y). Here, Q S (Y) oc e wc * y is the saturation momentum in the target, and it is also 
the transverse momentum scale which marks the transition from weak to strong scattering, for 
both the dipole and the quadrupole. 

The discussion of the dilute regime is straightforward. Introducing the dipole T-matrix 
operator T XlX2 = 1 — S X1X2 and the corresponding scattering amplitude (T x1X2 )y, then for 
r = \xi — xi\ <C l/Qs the scattering is weak, (T XiX2 )y "C 1, and the amplitude obeys the 
linearized (in (T)y) version of Eq. (2.9), that is, the BFKL equation: 



d(T XlX2 ) Y _ Q_ 
dY ~ 2vr 



■^ // ^xiX2z{T Xl z ~\~ T ZX2 T XiX2 )y • (3'1) 



Consider similarly the quadrupole: when r^j <C 1/Q a f° r an the six pairs (i, j) (with i, j = 
1,2,3,4) of external points, the scattering is necessarily weak, because there is no net color 
charge in the projectile over an area ~ l/Q 2 (which is the area where high gluon density effects 
become important in the target). Then, for consistency, it is enough to keep the lowest order 
terms, of order g 2 a 2 , in the perturbative expansion of 1— (Q)y- These are obtained by expanding 
the Wilson lines in Eq. (2.6) up to quadratic order. After also comparing with the corresponding 
expansion of S in Eq. (2.5), one finds 

1 (Q x\X2X2,x$)y — {T XlX2 T XlX[i + T XlX4 + T X2X3 T X2Xi -f- T x ^ X4 ^y j (3-2) 

where it is understood that, in the r.h.s., (T)y is the dipole amplitude in the dilute regime and 
obeys the BFKL equation (3.1). Eq. (3.2) also tells us that, for generic configurations at least, 
the quadrupole scattering can become strong when at least one (which necessarily means at 
least three) of the six transverse distances is of order 1/Q S , or larger. However, care must 
be taken when discussing very asymmetric configurations (see below). 
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Consider now the strong scattering regime at 3> 1/Q S (Y), where we shall study the 
approach towards the 'black disk' limit, i.e. the way how ^-matrices like (S)y and (Q)y 
approach to zero. The corresponding analysis for the dipole is well known [35, 36], but it will 
be briefly reviewed here, in preparation for the quadrupole and higher n-point functions. The 
crucial observation is that, at least for large N c , this approach is controlled by the virtual terms 
in the respective evolution equations 2 . Indeed, at large N c we have, e.g., (SS)y — (S)y(S)y, and 
then the 'real' term in the r.h.s. of Eq. (2.9), which is quadratic in (S)y, is much smaller than the 
linear, 'virtual', term there whenever (S)y <C 1. More precisely, when r = \x\ — X2\ S> 1/Q S (Y), 
the linear term dominates over the quadratic one unless either \x\ — z\ or \z — x 2 \ is as small 
as 1/Q S (Y). Then, the r.h.s. of Eq. (2.9) is dominated by the logarithmic regions of integration 
where one has 

1/Q B (Y) < \x! -z\, \z-x 2 \ «r. (3.3) 

This discussion implies that the approach of the dipole S-matrix towards the 'black disk' limit 
is determined by the following, simplified, equation 



d\n(S(r)) Y __ f d£ 2 
8Y - a Jun2 z 2 



a I ri_ = -aln[r 2 Q 2 s (Y)], (3.4) 



which is correct to leading logarithmic accuracy with respect to the large logarithm ln(r 2 Q 2 ). 
That is, in writing the r.h.s. of Eq. (3.4), we control the coefficient in front of the logarithm but 
not also the constant term under the logarithm. Using the fixed coupling relation ln[r 2 Q 2 (^)] — 
uia(Y — Yq), where Yq is such that Q s (Yq) = 1/r, we finally arrive at 

(S(r)) Y * (S(r)) Yo expj - i-ln 2 [r 2 Q 2 (y)]|. (3.5) 

This expression is valid for Y > Yq, where Yq is roughly the rapidity at which the scattering of 
the dipole with size r becomes strong. 

We now turn to the quadrupole and similarly study the strong scattering regime for generic 
configurations of the four external points in the transverse plane. Note however that our ultimate 
objective in this study is different, and actually more ambitious, than in the corresponding study 
of the dipole: what we are truly interested in, is not the law for the approach towards the 
black disk limit (although this law will emerge too from our subsequent analysis), but rather 
the functional relation between the quadrupole 5-matrix (Q)y an d the dipole one (S)y, as 
predicted by the solutions to the respective evolution equations with the virtual terms alone. 
Indeed, as announced in the Introduction and will be explicitly checked on the final results, this 
functional form is correct also outside the strong scattering regime — namely, it has the right 
limit at weak scattering, as shown in Eq. (3.2). 

Note first that a quadrupole configuration is bound to be in the strong scattering regime 
whenever all the four quark-antiquark separations (r\ 2 , ?"i4, ?"23 and ^4) are much larger than 



2 This statement is strictly correct only within the limits of the JIMWLK evolution. As emphasized in [36], 
the inclusion of Pomeron loops [37-40] in the analysis would modify the approach towards the black disk limit, 
even for large N c . Still, the functional form of this approach, as shown in Eq. (3.5), should remain correct: only 
the overall coefficient in front of the logarithm squared should be reduced by a factor of 2 [36] . 
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X3 x 4 



(a) 



(b) 



(c) 



X 2 X4 X3 

(d) 



Figure 1. Special quadrupolc configurations that we shall use to illustrate our general results. The dots 
represent the positions of the fermions in the transverse plane. (A bar on a transverse coordinate refers 
to an antiquark.) The lines have no physical meaning, they are meant only to delimitate the shape of 
the figures, in conformity with the names that we use to refer to them. 



1/Qs, independently of the relative separations between the two quarks (ris) and the two an- 
tiquarks (724). Indeed, when r q q 3> 1/Q S , there is uncompensated color charge over distances 
> 1/Qs, which implies that the scattering is strong. In this regime and for large 7V C , the evo- 
lution towards the unitarity limit is controlled by the virtual terms, by the same arguments 
as for the dipole. After neglecting the 'real' terms, Eq. (2.10) simplifies considerably and, in 
particular, it becomes 'diagonal' with respect to the quadrupole configuration: both the l.h.s. 
and the r.h.s. involve the same configuration {x\, X2, X3, x^). This makes it possible to study 
the approach towards the unitarity limit configuration by configuration and thus obtain explicit, 
analytic solutions. 

The case of a general configuration (which respects the conditions stated above) will be 
addressed in the next section. But before that general discussion, it is instructive to consider 
a simple, particular case as a warm-up. In what follows, we shall often use the four special 
configurations shown in Fig. 1 in order to illustrate our results. These configurations are useful 
in that they involve only few independent transverse separations rjj, due to their high degree 
of symmetry. Besides, two of them — the 'line' in Fig. l.a and the first 'square' in Fig. l.b - 
have been previously chosen in the numerical studies in [22], so for them we know already the 
respective numerical predictions of the JIMWLK equation. This information will be useful for 
testing our forthcoming analytic results. 

We shall pick the 'line' configuration, where we place the two quarks on top of each other 
and the same for the antiquarks (see Fig. l.a), as our warm-up example. For this example, we 
shall perform in detail all the steps to be repeated for the general case in the next section. (The 
three other configurations in Fig. 1 will be studied in Sect. 4, as special limits of our general 
result for the quadrupole.) As obvious from Fig. l.a, the 'line' configuration involves a single 
transverse scale r, which is the common size of the four qq pairs: r±2 = = r23 = = r. 
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When applied to this configuration, the general evolution equation (2.10) reduces to 
d(Q 

X1X2X1X2/Y ^* 



dY 



/ J^x\X2Z \{,S XiZ )y {Qzx-2XiX2)y (S ZX2 )y (,Qx\ZX\X2 ) Y 
K Jz 

-(Q x\x2x\x2)y )y] ■ 

(3.6) 



When the separation r is much larger than 1/Q S , the virtual terms (the last two terms in 
Eq. (3.6)) dominate and give 

c -2a Hr^Ql) [<Q(r))y + (S(r)) Y ] . (3.7) 

We have written the argument of the quadrupole simply as Q(r), because r is the only indepen- 
dent transverse scale for this particular configuration. After also using Eq. (3.4), it turns out 
that this inhomogeneous equation admits a relatively simple solution, which is (for Y > Yq) 



(S(r)) Yo (S(r))fa 



('line' in Fig. La) . (3.8) 



In view of its derivation above, one may expect this formula to hold only deeply at saturation, 
that is, for Y > Yq, with Yq the rapidity at which saturation sets in on the given transverse scale 
r. But as a matter of facts, Eq. (3.8) is more general than that: it has also the right limit at 
weak scattering, in agreement with Eq. (3.2), and hence it provides a global approximation. To 
see that, let us study the weak scattering limit of Eq. (3.8), by writing (S(r)) Y = 1 — (T(r))y, 
with (T(r))y <C 1, and then expanding the r.h.s. to linear order in (T(r))y. After doing that, 
one finds that the following relation is satisfied at Y, provided it was already satisfied at Yq : 

(Q(r))y ~ 1 - 4(T(r))y . (3.9) 

As anticipated, Eq. (3.9) is the correct weak-scattering limit for this particular configuration, 
as predicted by Eq. (3.2). 

Eq. (3.8) can be further simplified if one assumes that the initial conditions at Yq are 
provided by the MV model at large N c . Indeed, it is easy to check that Eq. (3.8) reduces to 

(Q(r)) Y = (S(r)) Y [l + In (S(r)) Y ] , (3.10) 

provided the above formula already holds at the initial rapidity Yq — a property which is indeed 
satisfied within the MV model at large -/V c . 

To summarize, by using Eq. (3.10) together with a reasonable approximation for the dipole 
S-matrix {S(r)) Y — say, the solution to the BK equation with initial conditions provided by 
the large- N c version of the MV model — one obtains a unified description for (i) the initial 
conditions at low energy, (ii) the BFKL regime at high energy, and (iii) the approach towards 
the black disk limit at all the energies. 

As demonstrated by the numerical analysis in Ref. [22], Eq. (3.10) agrees well with the 
numerical solution to the JIMWLK equation for this particular configuration. In the forthcom- 
ing section and in Appendix B, we shall generalize Eq. (3.8) to generic configurations for the 
quadrupole and the sextupole, and we shall also study some other particular configurations. 
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4 A global approximation for the quadrupole 



We shall now extend the analysis in the previous section to a generic configuration for the 
quadrupole — that is, we shall study the approximation obtained by keeping only the virtual 
terms in Eq. (2.10). In line with our general strategy, we shall rely on this approximation, 
which is strictly speaking valid only in the vicinity of the black disk limit, in order to derive 
a functional relation between the quadrupole and the dipole, that will be then promoted to a 
global approximation. 

After neglecting the 'real' terms and restricting the integrations over z to the regions yielding 
large logarithmic contributions, like in Eq. (3.3), we get 



d(Q Xl 



X2X3X4/Y 



dY 



a 

2 

a 

2 

a 

2 



[\n(rf 2 Q 2 s ) + ln(rteQs) + H^Ql) + H4 3 Q 2 S )] {Qx lX2 x a x 4 )Y 
[Hr 2 2 Q 2 ) + H4 A Ql) - ln(r 2 3 Q 2 ) - ln(r 2 4 Q 2 )] (S XlX2 ) Y (S x . 3X4 ) Y 

[ln(r 2 4 Q 2 ) + ln(r 2 2 3 Q 2 ) _ ln(r 2 3 Q 2 ) - ln(r 2 4 Q 2 )] (S x . 3X2 ) Y (S XlX4 ) Y , 

(4.1) 



where in writing the r.h.s. we have also factorized the 2-dipole S*-matrix, as appropriate at 
large N c . Now we can use Eq. (3.4) in order to express the logarithms in terms of the dipole 
S'-matrix. The merit of that is three-fold. First, we can express the quadrupole in terms of the 
dipole without worrying about the explicit solution for the latter. Second, in order to get the 
dominant logarithmic contribution we do the same kind of approximation in both the dipole 
and the quadrupole equations since their structures are the same. Therefore we probably have 
better control than originally announced, that is, we have better accuracy than the one given by 
the leading logarithmic terms. Third, we do not really need to specify the energy dependence 
of Q s and thus the approach is fully valid even if we consider a running coupling scenario (at 
least in the case where a is evaluated at Q s ). Thus, writing all logarithms in Eq. (4.1) terms of 
S we have 



9 {Qx\X2 X3X4 )y 

dY 



d 

^ln{S XlX2 ) Y (S XsX2 } Y (S XsX4 } Y (S XlX4 ) Y 



+ 



1 

+ 2 



d ^ (S XlX2 ) Y (S X3Xi ) Y 
®Y (S XlX3 ) Y (S X2X4 ) Y 

d , (S XlX4 ) Y (S X3X2 ) Y 
dY 



\SxiX3/Y (S: 



X2X4/Y 



(QxiX2XzX4)Y 

1X2IY {S XZX4 ) Y 
(S x \X4)y (S X3X2 ) Y . 



(4.2) 



This is an ordinary first order inhomogeneous differential equation whose general solution is 
straightforward to find; it reads 



{Qxi 



X2X3X4/ 



y = V (S XlX2 ) Y (S XsX2 ) Y (S X:jX4 ) Y (S XlX4 ) Y 



{Qxi 



X2X3X4/Y0 



-Y 



+ 



'Y 



dy 



(S Xi x2)y (S X3X2 ) Yo (S X3X4 } Yo (S XiX4 } Yo 

{S XlX3 }y{S X2X4 }y 8 (S XlX2 )y(S X3X4 )y + (S Xl x 4 ) y (S X3 x 2 ) y 



(S XlX2 )y(S X3X2 )y (S x 1x4)1 



dy 



(SxiX3)y{Sx2X4)i 



(4.3) 
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As explained in the Introduction, the above expression is guaranteed to be correct also for weak 
scattering, even though it has been obtained here via manipulations which are well justified only 
deeply at saturation. This is so since, at weak scattering, the relation between the quadrupole 
and the dipole becomes linear, cf. Eq. (3.2). Such a linear relation is then preserved by any 
evolution equation generated by a Hamiltonian — in particular, by the equations including only 
the virtual terms, which are generated, as we have seen, by a truncated version of the JIMWLK 
Hamiltonian (the first two terms in Eq. (2.2)). This argument about linearity is simple, but at 
the same time important for the present construction, so let us make a small digression in order 
to explain it better. 

Consider a set of three operators, A, B, and C, which obey a linear relation: A = C1B + C2C. 
Also assume that the evolution of these operators with 'time' Y is governed by some Hamiltonian 
H; e.g. dyA = HA. Then one can successively write 



was satisfied at the initial 'time' Yq. Note that specific form of H was irrelevant for the above 
argument; all that matters is the fact that this is a linear operator. 

Returning to the physics problem at hand, the relation (4.3) between the quadrupole and 
the dipole has been obtained by solving a Hamiltonian evolution equation. The corresponding 
Hamiltonian is incorrect in the dilute regime, where it would predict a wrong evolution for the 
quadrupole, or the dipole, taken separately. However, as a linear operator, it must preserve the 
linear relation (3.2) between the two operators Q and S which is valid in that regime. And 
indeed, it can be easily checked — by rewriting (S XiXj ) y = 1 — (T XtXj ) y and then expanding the 
r.h.s. of Eq. (4.3) to linear order in the various {T XiXj ) y s — that Eq. (4.3) predicts the correct 
relation, Eq. (3.2), between the quadrupole and the dipole at weak scattering and at rapidity Y 
provided this relation was already satisfied at the initial rapidity Yq. 

In view of the above, we conclude that Eq. (4.3), when used with a numerical solution 
of the BK equation (or some good approximation to it), provides a reliable approximation 
to the solution to the quadrupole equation which is valid in all regimes (at least for generic 
configurations, which are not very asymmetric). This solution provides a smooth (infinitely 
differentiable) interpolation between the BFKL solution in the weak scattering regime at rjj <C 
1/Qs and the correct approach towards the black disk limit at rjj ^> 1/Q S - 

The expression (4.3) has some other good properties. It is symmetric under the interchange 
of the two quarks (or the two antiquarks) at any Y, provided that this is true at the initial 
rapidity Yq. This is a true property of the JIMWLK evolution at large- N c , as it can be directly 
checked at the level of the evolution equations in Sect. 2. Furthermore, there is an interesting 
relation between Eq. (4.3) and the corresponding expression in the MV model [14, 20]. To see 
this connection, it is useful to write the dipole iv-matrix in the form 



Now, (4.3) becomes formally identical with the quadrupole formula in the MV model (for large 
N c ) provided one assumes Fy(xi,Xj) to be a separable function of Y and the transverse coor- 
dinates, plus an arbitrary function of Y. This property is certainly not satisfied in general, but 




(4.4) 



This equation implies that the linear relation A = c\B + c^C is preserved at any Y provided it 




(4.5) 



it is locally satisfied under some approximations — e.g., in the window for extended geometric 
scaling, where Ty(r) — ir 2 Q 2 {Y))^ s with ^y s ~ 0.63 (modulo more slowly varying logarithmic 
terms) [41-43], and deep at saturation where T Y (r) ~ {1 /2u) In 2 [r 2 Q 2 (Y)] (cf. Eq. (3.5)). This 
is also fulfilled in some dipole models, like the GBW model [44, 45], where IY(r) = r 2 Q 2 (Y)/A. 
Assuming this to be the case, then it is possible to explicitly perform the integral over y in 
Eq. (4.3) and thus obtain, after some algebraic manipulations, 



(Q 



XHB2X3X4/Y 



T Y (xi, x 2 ) + T Y (x 3 , x 4 ) - r Y (xi,xa) - T Y (x 2 , a? 4 ) 



+ 



r Y (xi,x 2 ) +r y (a; 3 ,a;4) 
ry(xi,a5 4 ) +T y (a;2,a;3) 



r Y (xi,x 4 ) - T Y {x 2 ,x 3 ) 
T Y (xi,x 3 ) - T Y (x 2 ,x 4 ) 



{S XiX2 )y (S X3X4 ) Y 

(S XlX4 ) Y (S X3X2 ) Y ,(A.6) 



rV(£Ci,£c 4 ) + T Y (x 2 ,x 3 ) - Y Y {x 1 ,x 2 ) - T Y (x 3 ,x A ) 

which is indeed the same as the corresponding expression in the MV model at large N c [14, 20] . 
Notice that the quadrupole at Yo does not appear anymore since in writing Eq. (4.6) we have 
assumed that the latter is already valid at Yq, as correct for initial conditions of the MV type. 

Returning to the general formula (4.3), let us notice that there are special quadrupole 
configurations, like those shown in Fig. 1, which have enough symmetry for the integral over 
y to be exactly doable, without additional assumptions. For instance, consider the case where 
all transverse separations between a quark and an antiquark are the same, that is r\ 2 = r 3A = 
r i4 = r 32 = r - (This includes the 'line', Fig. l.a, and also the first 'square', Fig. l.b.) Then the 
integrand in Eq. (4.3) is a total derivative, so one easily finds 



(Q)y 



(S(r)) Y 



(Q) 



Yo 



(S(r)) 2 



In 



(S(r)) 2 Y 



Y 



(S(r 13 )) Y (S(r 2i )) 



In 



0(r)) 2 



Yo 



Y 



(S(r 13 )) Yo (S{r 2i )) Yo 



(4.7) 



The arguments of {Q)y, which here is a function of r, r± 3 , and r 2 4, are kept implicit. 

The 'line' configuration in Fig. l.a corresponds to r\ 3 = r 2A = 0, and then Eq. (4.7) reduces 
to Eq. (3.10), as expected. As for the 'square' configuration in Fig. l.b, one has ri 3 = r 24 = V2r, 
so Eq. (4.7) gives 



(Q(r))y = (S{r)) 



(Q(r))Y 
(S(r)) 2 



+ In 



(S(r)) 2 



Y 



Yo 



(S(V2r)) 2 Y 



In 



(S(r)) 2 Yo 



(S(V2r)) 2 Yo 



Notice that the above reduce to 



(Q(r)) 



Y 



(S(r)) 2 Y 



1 + ln 



{S(r)) 2 Y 
(S(V2r)) 



Y 



('square' in Fig. l.b) . (4.8) 



(4.9) 



provided Eq. (4.9) is also valid at Yo, as it is indeed the case in the MV model at large- N c . More 
generally, Eq. (4.9) is valid at Y even if it was not valid at Yq provided the scattering of a dipole 
with the given size r was weak at Yo. (But of course, the scattering can be strong at Y > Yo for 
the same size r.) Indeed, assuming that (T(r))y ^ 1, and therefore also (T(^/2r)) Yo ^ 1, it is 
straightforward to show that the pieces linear in (T)y coming from the first and the last term 
of Eq. (4.8) mutually cancel and then Eq. (4.9) is again obtained. 

In Ref. [22], Eq. (4.9) has been compared to the numerical solution to JIMWLK equation 
and the agreement appears to be remarkably good — actually, even better than for the 'line' 
configuration. Clearly, expressions like (4.9) or (3.10), which involve a single transverse scale r, 
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exhibit geometric scaling in the same regime where the dipole does so. This feature too is in 
agreement with the numerical results in Ref. [22] 

A second class of configurations for which the integrand in Eq. (4.3) is a total derivative are 
those which are constrained by 7*13 = ri4 and r23 = r24, whereas the two other separations r\2 
and r34 are left arbitrary. The two last configurations in Fig. 1 — the second 'square', Fig. l.c, 
and the 'double triangle', Fig. l.d — belong to this category. For any configuration in this class, 
Eq. (4.3) yields a very simple, factorized, expression 



(those which are constrained). As before, Eq. (4.10) reduces to an even simpler expression 



provided this last formula holds already at Yq (a property which, once again, is satisfied within 
the MV model at large N c ). In particular, for the second 'square' configuration in Fig. l.c, one 
has ri3 = ri4 = r23 = r24 = r and r±2 = = \/2r, and Eq. (4.11) gives 



But the most interesting configuration in our opinion is the last one in Fig. 1: the 'double 
triangle'. What is remarkable about this configuration is the fact that the separation between 
the qq pair {x\, X2) and the second qq pair (2:3, 2:4) can be made arbitrary, yet the corresponding 
(S-matrix in Eq. (4.10) is independent of this separation. In particular, so long as r\% and r34 are 
much smaller than 1/Q S , the scattering of the quadrupole as a whole remains weak independently 
of the distance ~ ru separating the two pairs — including in the case where ru S> 1/Q S (Y). 
This may look counterintuitive since normally one expects the scattering to be strong whenever 
the separation between color charges is of order 1/Q S or larger. However, we believe that this 
result can be physically understood as follows: when the two quark-antiquark pairs (2:1,3:2) 
and respectively (2:3,0:4) are individually small, whereas the distance ~ ri4 between them is 
of order 1/Q S or larger, the color exchanges between the two pairs are strongly suppressed by 
saturation; then, the color charges compensates locally, within each of the small qq pairs, which 
therefore behave like two individual dipoles. In view of this argument, the factorized structure of 
Eq. (4.10) looks quite natural. It would be very interesting to verify this formula via numerical 
simulations of the JIMWLK evolution, along the lines in Ref. [22]. 
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A The evolution equation for the quadrupole 

In this section, we shall present a particularly streamlined derivation of the evolution equation 
for the quadrupole, using the 'dipole' version of the JIMWLK Hamiltonian in Eq. (2.2). (An 
alternative derivation can be found in Ref. [21].) 




(Q) Y = (S(r 12 ))Y(S(r 3 4))Y 



(4.11) 



(Q(t))y 



(S(V2r)) 2 Y 



('square' in Fig. l.c) . 



(4.12) 
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For pedagogy, we shall first derive the evolution equation for a color dipole, that is, the 
2-point function of the Wilson lines shown in Eq. (2.5). The action of the functional derivative 
on these Wilson lines is computed according to Eq. (2.8), thus yielding 

6^ J^fi v ^ Vx2 f 5 *^ 5 ™ tay k v ^ h + g 2 $ xl v5 X2U t b vlv X2 t a , (A.l) 

where in writing the r.h.s. we have dropped some terms proportional to d uv since the kernel 
■Muvz vanishes when u = v. 

Consider the first term in the parenthesis of Eq. (2.2), that is, the identity matrix. For it, 
both terms arising from Eq. (A.l) contribute equally and together yield the following contribu- 
tion to HS XlX2 : 



N£-l a 
~2Nf~ 2^ 



J^x\X2z S xi x 2 , (A. 2) 



where we have also used t a t a = — 1)/2N C . Next consider the contribution coming from the 
second term in the parenthesis of Eq. (2.2). After expressing adjoint Wilson lines in terms of 
fundamental ones according to 

it becomes straightforward to show that this second term yields the same contribution as the 
first one, that is, Eq. (A. 2). The third term in the Hamiltonian leads to contributions like 

{Vl) aC V* tv(t a 4 V X2 t b ) = tv{Vl t c V X2 V\ t c V z ) . (A.4) 

By using the Fierz identity 

ti(t a At a B) = -trAtrB- — tr(AB), (A.5) 

we find that the corresponding contribution to HS XlX2 , which also is equal to the one coming 
from the fourth term in Eq. (2.2), is 

J^ X \ X2 z I S XlZ S ZX2 Q S XlX2 I . (A. 6) 



2 2n l xi X2 z y^ X \z^z X2 jy 2 

z 

Putting all contributions together and averaging over the color field of the target, we arrive at 
the dipole equation shown in Eq. (2.9). 

Now let us turn our attention to the object of interest, the quadrupole operator, as defined 
in Eq. (2.6). The two functional derivatives can now act on 6 pairs of Wilson lines while the 
remaining two Wilson lines are simply 'spectators'. Just for the sake of illustration we shall give 
here some intermediate steps for the terms that arise when the 'active' pair is the one composed 
of the two quarks at x\ and 2:3. Acting on the respective Wilson lines with the first term of the 
Hamiltonian, we find 



9 2 



Stt 3 N c Jz 



M XlX3 Mt a vlv X2 t a vlv X4 ), 



-15- 



which after simple manipulations using the Fierz identity in Eq. (A. 5) leads to 

Ast Xl X3Z ( S Xl x 2 Sx^X 4 r, Qx\X2X%X 4 ) • (A.7) 



^ J " 'xix^z y^xiX2^xsx 4 



c 



Acting with the second term of the Hamiltonian we obtain 
j M Xl x zZ [{V^ c V^tr{t^ 



and by also making use of Eq. (A. 3) we are lead to 

■fi^XlXsZ [ Sx^X2 Sx\X 4 AT 2 Qx\X2X^X 4 ) • (A.8) 



^ J " l x 4 x A z y^x z x 2 ^x 1 x 4 

The action of the third term of the Hamiltonian gives 



<r 



16vr3Ar c j 



M XlXsZ (Vl) ac Vf tiit a VlVx 2 t b VlV X4 ) + (Vl) ac VZ b tr(t b VlVx2t a VlV X4/ 

■ J z 

leading to 

Finally the action of the fourth term in the Hamiltonian gives 

M X1X3Z \{V}) ac v£tr(t a VlVx2t b VlV X4 ) + {v}) ac V° b Mt b <V X2 t a VlV Xi )~ , 



2 

9 



16n 3 N c 

yielding the same result as shown in Eq. (A. 9). Putting together all the contributions, we see 
that the terms explicitly suppressed by cancel each other, as anticipated. Thus, the action 
of the Hamiltonian on the pair made with the two quarks at x\ and x% gives 

— j -^xix:jZ ^S Z X2 Qxizx^Xi S ' zx 4 Qx\X2X-3,z S X4 x2^x;iX 4 S X ^X2 Sx\x 4 ) • (A. 10) 

Similarly, when acting on the quark and the antiquark at x\ and X2 respectively, we obtain 

~. / J^x\X2Z \Sx\zQ zx2Xzx 4 ~\~ S Z X2 Qx\zxi,x 4 Qx\X2X3,x 4 S X4 X2 Sx%x 4 ) ■ (A.ll) 
4?r Jz 

The action of the Hamiltonian on the other possible pairs of Wilson lines is obtained via appro- 
priate permutations in either Eq. (A. 10) or Eq. (A.ll). Putting all terms together and averaging 
over the target field, we arrive at the evolution equation shown in Eq. (2.10) in the main text. 



B The sextupole 

In order to demonstrate the generality of our procedure, we shall succinctly present here the 
corresponding analysis for the next non-trivial high-point correlator, which is the sextupole 
defined in Eq. (2.7). The respective evolution equation has not been yet presented in the 
literature, thus we have derived it on this occasion. We shall not give here the details of the 
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derivation, since this is merely a lengthy but straightforward repetition of the manipulations 
shown in the previous Appendix. This is the final result: 



d? ~ 4^ 
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X2X4Z) {Q X5X6X1 zQ ZX2X3X4} Y 
X4X5Z 
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~\~ -M^xsxqz ~i~ ■A / tx(jXiz ) {^ XiX 2 X ^ X 4 X ^ X q }y 



X3X5Z I X*- 5 X3X4 ^XlX2X5Xfi/Y 
X4X§z) (S X5X ^ QxiX2X3Xq)y 
X4X§z) (Sx^xq Qx\X2X3X4)Y 

X2X6z){£>XiX(jQx3X4X$X2 )y\- (B.l) 

We notice that all the subleading terms of relative order l/AT? have canceled in the final equation, 
as expected for a single trace operator. As a cross check, one can see that when choosing, for 
example, x$ = xq, the above equation reduces to the quadrupole equation (2.10). 

The subsequent manipulations follow the general procedure outlined in Sect. 4. First, one 
separates the 'virtual' terms in the r.h.s. of Eq. (B.l), that is, the terms in which the integration 
variable z appears only in the dipole kernel, but not in the correlation functions; these are the 
terms proportional to (5"}y and to (SQ)y — (S)y(Q)y- Then in evaluating the integrals 
over z, one keeps only the dominant logarithmic contributions and use Eq. (3.4) to express 
the logarithms in term of the dipole derivative. We thus arrive at an ordinary first order 
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inhomogeneous differential equation whose solution is 



(6) 



X2X :i X4X 5 X 6 / 



y ~ V (Sx 1 X2)y(Sx3X2)y(Sx3X4}y(Sx 5 X4}y (S X5 xe)Y (S Xl X(i}Y 

4x2x3x4x5x5) Yq 

{S x 



+- 



Y 



Yo 



(S X ix 2 ) y (£>x 3 X2 )y (Sx s x 4 )y \S X5X4 )y (S X5 x e )y (S Xl xo )y 



id \ in \ —~\-n \^^i^)y\^x 3 xe)y 

\>JxiX2 1 y vxx 3 X4X^X(, }y n in- 



+ (Sx 
J r{S X 3X4)y 
J r{S X c,X4)y 
~\~{Sx 5X6 )y 



{Qx 1 X4X 5 x 6 )y 
{QxiX2X5Xe)y 



{Qxi 



X2X;iX6/y 



{Qxi 



X2XsX4/y 



I In 



J- in 



J- in 



J- in 



+ {Sxix^y (Q 



xe,/y\^x-sX4X 5 X2/y 



In 



(S X1X3 ) 


y{Sx2Xe) 


y 


(S XS X2 


*>y(S Xl x4 


>v 


(S Xl x 3 


)y{S X2 X4 


h 


(S XS X4 


)y{S X5 X2 




(Sx 2 X4 


"lyiSxsXs 


h 


{Sx 5 X4 


)y{S X3 x e 




(S X3 X5 


~>y(S X 4Xe 


h 


(Sx 5 x 6 


}y(S Xl x4 


>y 


(Sxixs 


)y(S X 4Xfi 


)y 


{S Xl x 6 


}y{S X5 X2 


)y 


{Sx±X5 


}y(S X 2X6 


)y 



(B.2) 



Again, one can check that choosing, for example, x§ = xq, the above solution reduces to the 
quadrupole one given in Eq. (4.3). Expanding Eq. (B.2) in the weak scattering region, we find 

1 {SxiX2X 3 X4X5X6 )y — {TxiX2 T Xl x 3 ~\~ T Xl X4 T Xl xs ~i~ T Xl xg T X2 X 3 T X2 X4 ~i~ -^S 



+ T X 



X 2 X 5 

^x 3 xs ~\~ T X3 xs TX4X5 ~ T XiX& + T X5X6 )y i (B-3) 



where we have also used Eq. (3.2) for (Q)y and we have assumed that Eq. (B.3) is already valid 
at Yq. It is straightforward to confirm that the same result is obtained when we expand Eq. (2.7) 
to order g 2 a 2 . 

As a particular example, we shall consider the 'line' configuration, already studied in the case 
of the quadrupole. This is obtained by putting all the quarks at the same point, x\ = x% = X5, 
and similarly for the anti-quarks, X2 = x^ = xq, with the two points separated by a distance 
r. (One can trivially visualize this configuration by looking at Fig. l.a and adding there 2:5 
and xq at the left and right ends, respectively.) Then by also using Eq. (3.8) we see that the 
y-integration in Eq. (B.2) can be exactly performed and gives 
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where clearly all quantities are evaluated at r. The above simplifies to 
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if we assume that the above and Eq. (3.10) are already valid at Yq or if we assume that the 
scattering at Yq for the given r is weak. The sextupole in Eq. (B.5) exhibits geometric scaling 
in the same regime where the dipole does so. 
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